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Abstract. We review recent developments in double-beta decay, focusing on 
what can be learned about the three light neutrinos in future experiments. We 
examine the effects of uncertainties in already measured neutrino parameters and 
in calculated nuclear matrix elements on the interpretation of upcoming double-beta 
decay measurements. We then review a number of proposed experiments. 
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1. Introduction 

Double-beta {(3(3) decay is a second-order weak process in which two neutrons inside 
a nucleus spontaneously transform into two protons. To conserve charge two electrons 
must be emitted. If lepton number is also to be conserved two antineutrinos must be 
emitted as well. This lepton-number-conserving process — /3/3(2u) decay — has been 
observed in several nuclei. Lepton number is not associated with a gauge symmetry, 
however, and so its conservation is not sacrosanct. If lepton number is violated, e.g. 
through the propagation of Majorana neutrinos, then a variant of the decay in which 
no neutrinos are emitted — (3f3{^v) decay — may also occur, though it has never been 
observed. In (3f3{^v) decay a virtual neutrino is emitted by one neutron and absorbed 
by the second. (3f3{^v) decay can also occur through the exchange of other particles, 
perhaps some predicted by supersymmetric models, between the nucleons. Although we 
address this and other more exotic possibilities here, we are most interested in (3(3{Qu) 
decay mediated by light Majorana neutrinos because its rate depends on the absolute 
neutrino mass scale, a number on which we currently have only a generous upper limit. 

In recent years experimenters have discovered that neutrinos have nonzero masses 
and mixings, and have pinned down many of the associated parameters. But although 
we now know the differences between the squares of the masses, we do not know the mass 
of the lightest neutrino, nor the pattern in which masses of the three active neutrinos 
are arranged. /3/?(0z/) experiments have the potential to teach us about these matters, 
and in this review we focus on the question of how much we can expect to learn from 
experiments in the next decade. 

Previous reviews have thoroughly addressed other aspects of (3(3 decay. Early papers 
by Primakoff and Rosen (1959), Haxton and Stephenson (1984), and Doi et at (1985), as 
well as more recent reviews by Tomoda (1991), Suhonen and Civitarese (1998), Vergados 
(2000) and Klapdor-Kleingrothaus (2000) presented the theoretical formalism in great 
detail. A comprehensive review by Faessler and Simkovic (1998) in this journal devoted 
particular attention to (3j3{^v) decay in supersymmetric models and to the calculation of 
the nuclear matrix elements governing all kinds of (3(3 decay in the Quasiparticle Random 
Phase Approximation (QRPA). We will treat these topics, particularly the first, more 
briefly. A recent review by Elliott and Vogel (2002) examined the impact of recent 
discoveries in neutrino physics on (3(3{Qv) decay, but since then several experiments (e.g. 
KamLAND and WMAP) have reported results and a number of papers interpreting 
the results have appeared, with the consequence that neutrino masses and mixings are 
better understood. We will focus most intently on the additional neutrino physics that 
can be extracted from (3(3(0u) decay in light of these very recent results. The discussion 
will require us to examine the accuracy with which the (3(3(Qv) nuclear matrix elements 
can be calculated, in addition to surveying present and future experiments. 
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2. General Theory of (3(3 Decay 

This subject has been covered extensively in the reviews listed above and we do not 
present it in great detail here. Since we focus on what can be learned about neutrinos 
from /9/9(0z/) decay, we must begin with a brief discussion of Majorana particles. 

We define the left and right handed components of a Dirac 4-spinor by ^l,r = 
[(1 =f 7 5 )/2]^. In the Standard Model, only left-handed neutrinos interact. Because 
neutrinos are neutral, however, there is an additional way to construct left-handed 
neutrino fields. The charge conjugate field \l/ c = i^y 2 \]/* is also neutral; neutrino fields 
can be linear combinations of ^ and \& c and, since {^r) c is left-handed (and {^l) c is 
right-handed) we can define independent Majorana neutrino fields that are their own 
charge conjugates (antiparticles) via 

*L + (*L) C y ^R+^rY m 

Mass terms in the Lagrangian can couple these two kinds of fields to themselves and 
each other; the most general mass term has the form 

C M = -M L vv - M R XX - M D (uX + Xv) . (2) 

The parameters Ml and Mr are "Majorana" masses for the v and X fields, and Mr, is 
a "Dirac" mass that couples the two. For iV flavors of neutrinos these masses can be 
arranged in a matrix: 



where the matrices M.l, M-r, and M.r> are now N x N. If M. L and M.r are zero, the 
Majorana z/'s pair with an X to form iV Dirac neutrinos. At another extreme, in the 
"see-saw" mechanism (Gell-Mann et al 1979; Yanagida, 1979; Mohapatra and Senjanovic 
1980) — here with one flavor, for simplicity — M L « M D « Mr and the resulting 
eigenstates are Majorana neutrinos, the lightest of which has mass m ~ M^/Mr. 

In general the eigenstates of M. represent Majorana neutrinos that are related to 
the fields v and X by the "mixing" matrix of eigenvectors: 

where U and V are iV x 2N submatrices. We assume here that the states $ are either 
absent or extremely heavy, in which case we can work with a nearly unitary N x N 
section of U that mixes the light neutrinos: 

V l - J2 U lm<Pm ■ (5) 
m 

The original states v\ with definite flavor that enter the Lagrangian are linear 
combinations of the states <p m with definite mass. 

The matrix U nominally has iV 2 parameters, iV(iV — l)/2 angles and iV(iV + l)/2 
phases. iV of the phases are unphysical (Kobazarev et al 1980), leaving iV(iV — l)/2 



CONTENTS 



5 



independent physical phases. For three active neutrinos, the mixing matrix can be 
written in the form 

( C12C13 S12C13 s 13 e' lS \ 



U 



~ S12C23 — Ci2S23Si3e J<5 C12C23 — Si2S23Si3e J<S S23 c 13 

\ S12S23 - c i2 c 23 s 13 e lS -C12S23 - si 2 C2 3 si 3 e l5 c 23 ci 3 j 



diag{e * ,e » ,1} , (6) 



where and stand for the sine and cosine of the angles dij, 5 is a "Dirac" phase 
analogous to the phase in the CKM matrix, and the other two phases ot\ and a 2 affect 
only Majorana particles. (They can't be rotated away because Majorana neutrinos don't 
conserve lepton number.) 

The neutrino masses and mixing matrix figure in the rate of neutrino-mediated 
/3/3(0v) decay. If the weak quark-lepton effective low-energy Lagrangian has the usual 
V-A form, then the rate for that process is 

[T^]- 1 = E / \ZoA 2 S(E e i + E e2 - A£)|^|^ , (7) 
spins 

where Z 0u is the amplitude for the process and AE is the Q-value of the decay. 

The amplitude Zq v is evaluated in second-order perturbation theory and can be 
written as a lepton part contracted with a hadron part. The lepton part of the amplitude, 
containing outgoing electrons and exchanged virtual Majorana neutrinos of mass rrij, 
emitted and absorbed with amplitude U e j, is 

- J /E (gl^-^Wl " Ta)^^ 1 " ">^°M V% , (8) 

where q is the 4-momentum transfer. The term with q p vanishes and the mj in the 
denominator can be neglected for light neutrinos, so that the amplitude is proportional 
to 



m x \U el \ 2 + m 2 \U e2 \ 2 e [{a2 - ai) + \U e3 \ 2 e^~ ai - 2S ^ .(9) 



The absolute value has been inserted for convenience, since the quantity inside it is 
squared in equation (JJj) and is complex if CP is violated. 

The hadronic part of the amplitude must be evaluated between initial and final 
nuclear ground states, with the states in the intermediate nucleus summed over in the 
same way the virtual neutrino's momentum must be integrated over. For the amplitude 
to be appreciable, the wavelength of the virtual neutrino cannot be more than a few 
times larger than the nuclear radius R, i.e. only momenta q £ 1/R ~ 100 MeV will 
contribute. The excitation energies of the intermediate nuclear states generated by the 
hadronic part of the amplitude are all significantly smaller and so to good approximation 
the individual energies of those states in the "energy denominator" can be neglected or 
replaced by an average value E (to which the expression is not very sensitive) and the 
states can be summed over in closure. The result, in the allowed approximation for the 
weak hadronic current, is 
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,2 



M GT _9V M F 

9a 



2 



2 



Go^A^^lMo.Hm^) 2 . (10) 



Here Gq v (AE, Z) comes from the phase-space integral (which depends on the 
nuclear charge Z through the wave functions of the outgoing electrons), qa and gy 
are the weak axial-vector and vector coupling constants, and the M's, which are nuclear 
matrix elements of Gamow- Teller-like and Fermi-like two-body operators, are defined as 

Ml =(f\E H (r jk ,E)T+T+\i) (11) 

Mg T = {f\Y,H{r jk ,E)a r a k T?T+\i) . (12) 

The function H, which depends on the distance between nucleons and (quite weakly) on 

the average nuclear excitation energy in the intermediate nucleus, is sometimes called a 

"neutrino potential" and has approximate form 

. -. 2R f°° , q sin qr , . 

H(r, E) = — / dq— - -. H — - , 13 

V ' ; vrr Jo H u(u + E-[M i + M f }/2) ' K ' 

where Mi and Mf are the masses of the initial and final nuclei. 

The approximate expression equation (fTU|) is typically accurate to within about 
30%, the largest corrections coming from "induced" terms (weak magnetism, induced 
pseudoscalar) in the hadronic current (Simkovic et al 1999). For /3(3(2v) decay, the 
expression for the rate [T^] -1 is similar to equation (jlOj) . the differences being in 
the phase space factor G^^AE 1 , Z) ^ Gq u (AE, Z) and the matrix elements M<[ u and 
M2J ', which don't contain the neutrino potential H but do contain energy denominators 
because, without the intermediate neutrino, closure is not a good approximation. For 
later reference, we give the relevant (3f3{2v) matrix elements here: 

mi = v (f\^-rnn^^) (14) 

2u E m -[Mi + M f ]/2 K } 

m gt _ xr (f\Zj ff j T ?\ m )( m \^kVkT£\i) 

2 » ~\t E m - [Mi + M f )/2 ■ 

Other mechanisms besides light-neutrino exchange can drive 0[3(Ou) decay, and 
we discuss them briefly later. No matter what the exchanged particles, however, the 
occurrence of f3(3(0i>) decay implies that neutrinos are Majorana particles with nonzero 
mass (Schechter and Valle 1982). 



3. Phenomenology of Neutrino Properties and Double-Beta Decay 

The remarkably successful worldwide neutrino physics program has revealed much about 
neutrinos over the past decade. We now know that they mix and we have initial values 
for the mixing matrix elements. We know the differences between the squares of the 
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neutrino masses and the number of light active neutrino species. There is much we 
still don't know, however. In this section we summarize what we have learned and its 
implications for /3/3decay experiments that seek to learn more, focusing as mentioned 
above on the exchange of three light species of neutrinos. Other (5(5 possibilities are 
discussed later. 

What aspects of still-unknown neutrino physics is it most important to explore? 
Although the answer is to a certain degree a matter of opinion, it is clear that the 
absolute mass scale and whether the neutrino is a Majorana or Dirac particle are crucial 
issues. (5(5 decay is t he only laboratory process that can test the absolute mass scale with 
a sensitivity near \j <5m 2 tm (defined and discussed below). More importantly, whether 
the neutrino is Majorana or Dirac is a completely open question, and (5(5 decay is the 
only practical way to address it. Because future /3/3-decay experiments will be sensitive 
to a range of masses that includes \J 5m 2 tm and therefore at least one neutrino, even 
null results will have significant impact on our understanding. We believe that (5(5(0v) 
decay should be part of any future experimental neutrino program. 

3.1. Neutrino- Oscillation Parameters 

In this subsection, we define the various neutrino-oscillation parameters and discuss 
the ways their uncertainties affect our ability to extract (mpp) from a (5(5 experiment. 
The oscillation experiments have provided data on the mixing-matrix elements and 
the differences in the squares of the mass eigenvalues (Smfj = (m 2 — ml)). From the 
atmospheric-neutrino data, we have \Sml 3 \ = <5m 2 tm = 2.0tl'® x 10~ 3 eV 2 (90% CL) and 
# atm = ~ 45 degrees (Saji 2004). The combined results of solar-neutrino experiments 
and reactor experiments (Ahmed et al 2004) give 5m\ 2 = 5ml ol = 7.llJ;g x 10~ 5 eV 2 
(much less than <5m 2 tm ) and 6\2 = #soi = 32.5l 2 ; 3 degrees (68% CL). From reactor 
experiments, we have the limit #13 < 9 degrees (68% CL) (Hagiwara 2002). Note that 
other authors obtain modestly different results (see, e.g. Bahcall and Pena-Garay 2003, 
de Holanda and Smirnov 2003, Smy et al 2004) for the solar-reactor parameters. The 
Super-Kamiokande collaboration has also presented an independent value for |5m| 3 | of 
2.4 x 10~ 3 eV 2 (Ishitsuka 2004). Finally the limits on 13 depend on the specific value 
of 5ml 3 . For these reasons and the fact that the precise values for these parameters are 
rapidly changing, the discussion to follow should be considered illustrative only. 

The sign of <5m 2 ol is known; the lighter of the two eigenstates participating 
significantly in the solar oscillations, which we call v\, has the largest v e component. 
(The third eigenstate z/ 3 contains very little is e .) We know that v\ is slightly lighter than 
V21 but we don't know whether z/3 is heavier or lighter than this pair. If z/ 3 is heavier, the 
arrangement of masses is called the "normal hierarchy" (with two light neutrinos and 
a third significantly heavier one); if it's lighter the arrangement is called the "inverted 
hierarchy" . When all three masses are significantly larger than yj <5m 2 tm the hierarchy 
is referred to as "quasidegenerate" , no matter which eigenstate is the lightest. One of 
the largest questions left in the neutrino world is "Which of the hierarchies is realized 
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Minimum Neutrino Mass (meV) 

Figure 1. The effective Majorana mass (nipp) as a function of the lightest neutrino 
mass. 

in nature?". The answer will tell us about the overall scale of the neutrino masses as 
well as the order in which they are arranged. 

The central values of the oscillation parameters and Eqn. El determine a range of 
(m/3/3) values for a given value of mi. Many authors have analyzed the dependence (see, 
e.g., Vissani 1999, Bilenky et al 1999, Klapdor-Kleingrothaus et al 2001, Matsuda et 
al 2001, Czakon et al 2001, Elliott and Vogel 2002, Feruglio et al 2002, Giunti 2003, 
Joaquim 2003, Pascoli and Petcov 2003, 2004, Sugiyama 2003, Bilenky et al 2004, 
Murayama and Pena-Garay 2004) and Fig. Q shows the results. The bands indicate the 
range of possible values, which depend on the unknown phases in the mixing matrix. 
The borders indicate CP-conserving values of the phases, e^ ai ~ a ' 2 ^ = ±1. 

The observation of [3(3{Qv) decay would have profound implications regardless of 
uncertainty in the deduced value of (m^). It would show that neutrinos are massive 
Majorana particles and that the total lepton number is not a conserved quantity. But 
it is interesting to ask whether one can use a measurement of (m^) to discern the 
correct hierarchy. At high values of the minimum neutrino mass, the mass spectrum is 
quasi-degenerate, and the bands in Figure d are not resolved. For a minimum neutrino 
mass of about 50 meV, the degenerate band splits into two, representing the normal 
(mi lightest) and inverted (m^ lightest) hierarchies. Figure Q appears to imply that it 
would be straight-forward to identify the appropriate band at these low mass values. 
However, there are uncertainties in the oscillation parameters and the matrix elements 
that are not represented in the figure. 

One can address the question of how difficult it would be to distinguish the 
hierarchies by comparing the maximum value for the normal hierarchy ((m^)^™) 
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with the minimum value for the inverted hierarchy {{mpp) 1 ^^) that can result from 
the parameter uncertainties. When the lightest neutrino mass is small, the expressions 
for these values are simple (Pascoli and Petcov 2003). When mi is near zero, (mpp)^^ 
occurs for constructive interference between the contributions from the and terms 
in Eq. El when CP is conserved: 



From Eq. [H3 it is clear that (mpp)^^ is maximal when 9i 3 , 9 so \, and the <5m 2 's are as 
large as they can be. 

( m A3)mm i s minimal with the same conditions on #13 and 9 so \, but the smallest 
allowed value for 5m 2 t ■ 



If we use the appropriate extremum values for the oscillation parameters in Eqns. 
and EH we find (mpp)^ « (9.1 meV) (0.327) (0.976) + (55 meV)(0.024) = 4 meV 
and {mpp) 1 ^ ~ (36 meV) (0.345) (0.976) = 12 meV. These numbers are sufficiently 
different, at least when using our low-confidence-level uncertainty ranges, that one could 
discriminate between the two solutions. 

Since the precision of the oscillation parameters is likely to improve with future 
experiments, they would not appear to be a primary concern. Even so, it's interesting to 
see the effects of uncertainties in individual parameters. These effects can be determined 
by propagating the uncertainty in each parameter through to the uncertainty in (mpp). 
The results are shown in Tabled It is clear that #13 affects (mpp)^ x a great deal. It's 
also clear that Sml tm is critical for (mpp) 1 ™^. Finally, SO \ is important for both. In 
short, improved precision for 6*13, 9 so \ and Sm1 tm would help with the interpretation of 
a (3(3(Qv) experiment. 

The parameter 9 so \ has an effect in Figure^ that we haven't yet mentioned. If it 
has the right value, cancellation can drive {mpp) to very small values. But as long as 
solar mixing is substantially different from maximal, the cancellation is possible only 
over a narrow range of values for the lightest mass, and complete cancellation is not 
possible at all in the inverted hierarchy. 

The uncertainty in \M 0v \ has been a source of concern for a long time. Typically, 
it has been assumed to contribute of a factor of 2-3 times (mpp) to the uncertainty in 
(mpp). This uncertainty clearly dwarfs any from the oscillation parameters and thus is 
the primary issue. We address it later. 

3.2. The Absolute Mass Scale 

As already mentioned, we know that at least one neutrino has a mass greater than 



V ^ m atm ~ 45 meV; sensitivity to this value is therefore the goal of most future neutrino- 
mass experiments. In this subsection we briefly compare the potential of other mass 
measurements. 




(16) 
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Table 1. A summary of the impact on (mpp) in the normal and inverted hierarchies 
of the oscillation-parameter uncertainties. For the central values of the parameters, 
( m /3/3)max an d ( m /3p}min are 2-4 meV and 19 meV, respectively. See Pascoli and Petcov 



(2003) for a similar analysis. 


Oscillation 


Parameter 


Range 




Range 


Parameter 


Range 


in 


in 






8.1 - 9.1 meV 


2.3 - 2.6 meV 




N.A. 


\JSm 2 atm 


36 - 55 meV 


3.2 - 3.7 meV 


15.2 


- 23.2 meV 






(with #i3 = 9°) 






#sol 


30.1 - 34.9 cleg 


2.1 - 2.7 meV 


15.5 


- 22.4 meV 


6*13 


- 9 cleg 


2.4 - 3.5 meV 


18.6 


- 19.0 meV 



There are many ways to measure the mass of the neutrino, (see Bilenky et al 
2003) for a nice summary of the techniques). The best of these are /3/3(0z/)decay, (3 
decay, and cosmological observations. These three approaches are complementary in 
that they determine different combinations of the mass eigenvalues and mixing-matrix 
parameters (see Equations I23|) . A measurement of f3/3(0u) decay determines a coherent 
sum of the Majorana neutrino masses because {mpp) arises from exchange of a virtual 
neutrino. Beta decay measures an incoherent sum because a real neutrino is emitted. 
The cosmology experiments measure the density of neutrinos and thus a parameter 
proportional to the sum of the neutrino masses. 

The present limit (mp) < 2200 meV (95%CL) comes from tritium beta decay 
(Lobashev et al 1999, and Weinheimer et al 1999). This limit, when combined with the 
oscillation results, indicates that for at least one neutrino: 

45 meV <m 8 < 2200 meV . (18) 

The two /5-decay experimental groups just referenced have joined forces to form 
the KATRIN collaboration (Osipowicz et al 2001). They propose to build a large 
spectrometer that exploits the strengths of both previous efforts. KATRIN hopes to 
reach a sensitivity to (mp) near 200 meV (Bornschein 2003). 

Massive neutrinos would contribute to the cosmological matter density (Hannestad 
2003) an amount, 

tt u h 2 = S/92.5 eV , (19) 

where Q is the neutrino mass density relative to the critical density, lOO/i is the Hubble 
constant in km/s/Mpc, and S = mi + m 2 + is the sum of the neutrino masses. 
The neutrinos are light, however, and cluster with cold dark components of the matter 
density only for scales larger than 

k ~ 0.03y/(m u /leV)n m h Mpc -1 . (20) 

For smaller values of k, perturbations are suppressed, and so measurements of the large 
scale structure (LSS) can provide constraints on the neutrino mass. Such constraints 
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are rather weak, though, unless used in conjunction with precise determinations of the 
other various cosmological parameters, which also reflect the size of perturbations. Most 
recently the WMAP collaboration (Bennett et al 2003) has provided precise cosmological 
data that supplement LSS data from the 2dF galaxy survey (Elgar0y et al 2002), CBI 
(Person et al 2002), ACBAR (Kuo 2004) and the Lyman-a forest data (Croft 2002). 
These data have been used in various combinations to derive 95% CL limits on E of 
< 0.75 eV (Barger et al 2003), < 1.0 eV (Hannestad 2003a), < 1.7 eV (Tegmark et 
al 2003), < 0.69 eV (Spergel 2003), and < 1.0 eV (Crotty et al 2004). There is at 
least one claim for a nonzero value: £ = 0.64 eV (Allen 2003). An interesting paper 
by Elgar0y and Lahav (2003) points out the impact of the prior distributions on the 
resulting neutrino mass limit. 

Future measurements by the Sloan Digital Sky survey (SDSS 2003) and the 
PLANCK satellite (PLANCK 2003) may obtain limits on £ as low as 40 meV (Hu 1999). 
But the determination of E from cosmology is clearly complicated by the large number 
of correlated parameters that must be measured. Clean laboratory measurements of the 
neutrino mass will always be desirable. Ordinary /?-decay will be hard pressed to reach 
a sensitivity to \j5m1 tm , so (3(3(0u) experiments are especially important. 

3.3. The Majorana Phases 

Equation El shows the effective Majorana neutrino mass and its relation to the Majorana 
phases. When these relative phases are an integer multiple of 7r CP is conserved. In 
principle, the two relative phases (a = (a^ — ai), (3 = (— a.\ — 26)) have measurable 
consequences. In practice, determining them will be difficult. In this subsection, we 
discuss the physics of these phases. Many authors have examined the potential to 
combine measurements from (3(3 decay, tritium (3 decay, and cosmology to determine 
the Majorana phases. (See for example: Sugiyama 2003; Abada and Bhattacharyya 
2003; Pascoli et al 2002, Pascoli et al 2002a.) We can illustrate the main ideas through 
a simplified set of hypothetical measurements. Figure El shows a two-neutrino-species 
example of such a set. We took the mixing matrix and 5m 2 to be the best fit to the 
solar-neutrino data, with an arbitrary value for the Majorana phase a (of which there 
is only one) of 2.5 radians. We then made up values for E, (mpp), and (mp) assuming 
them to be the results of pretend measurements. Each curve in the m 2 vs. nil graph is 
defined by one of these measurements. We chose the value of E (from cosmology) to be 
600 meV, corresponding to a quasidegenerate hierarchy, and let (771/3) = 300 meV and 
( m 0p) — 171 meV. The rri2 versus mi curves from the "measurements" of the oscillation 
parameters, E, and (3 decay are: 

m 2 = S — mi 

m 2 = \Jm\ + b~m\ 2 

m 2 = ^{(m p )/U e2 y-{miU e i/U e2 y . (21) 
The (3(3 constraint is a little more complicated than that from (3 decay. The curve 
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is also an ellipse but rotated with respect to the axes. The constraint is given by the 
solution to the quadratic equation resulting from the two-neutrino version of Equation 

El 

a =Ut 2 

b = 2m 1 U 2 1 U 2 2 cos(a;) 
c = {rrixUlx) 1 - (m ee ) 

m 2 = (-b ± Vb 2 - 4ac)/2a (22) 

All of these equations express m 2 in terms of m\ and measured parameters and all should 
intersect at one (mi,m 2 ) point. However, because the point is overdetermined, the (3(3 
ellipse will intersect only for a correct choice of a. This provides a way to determine 
a. In Figure El we drew the (3(3 ellipse for a = 2.0 radians and for the "true" value of 
2.5 radians to show how the intersection does indeed depend on a correct choice of the 
phase. 

Although this two-species example is illuminating, it is overly simplistic. One needs 
to consider three species and two phases. In this case, the E constraint becomes a plane 
and the (3(3- and /3-decay constraints become ellipsoids. Figure El shows a model three- 
species analysis with E = 700 meV, 5m 32 positive, (mp) = 232 meV, and (mpp) = 159 
meV. (The phases were taken to be 2.0 and 2.5 radians, and U e s was taken to be = 
0.03.) The surfaces, shown in Figure El are defined by the equations for E, oscillations, 
(mp) and (mpp). Respectively, these are: 

E = mi + m 2 + m 3 

Sm 21 = m\ — m\ 

Sm 32 = ml — m\ 

(mp) 2 = m\U 2 el + m\Ul 2 + m 2 3 U 2 3 

(m ee ) 2 = mlU^ + m\Ul 2 + m\U^ + 2m l m 2 U 2 1 U 2 2 cos(a) 

+ 2m 1 m 3 U 2 1 U 2 3 cos(p) + 2m 2 m 3 U 2 2 U 2 3 cos(a + (3) . (23) 

These surfaces intersect at a point but (3(3 decay is the only measurement of those 
used that is sensitive to the phases. Thus a second pair of phases will also produce a 
consistent result. (That additional ellipsoid is not shown here.) Two experiments that 
depend on the phases are required to unambiguously determine both. Furthermore, to 
keep the plots legible, this analysis ignores any uncertainty in the measured parameters. 
The uncertainties will be large, at least initially, and therefore conclusions about the 
phases will be weakened. Barger et al (2002), for example, have noted that a significant 
improvement in \Mq u \ is necessary before the interpretation can be successful. The 
articles listed above consider the uncertainties in their analyses. 

When heavy right-handed Majorana neutrinos decay, they will violate lepton 
number. In the early universe, these decays would be out of equilibrium and could 
violate CP. The resulting net lepton number could be transferred to a net baryon 
number through standard weak interactions. Thus, this leptogenesis process (Fukugita 



CONTENTS 



13 



and Yanagida 1986) could explain the baryon asymmetry in the Universe. In principle, 
leptogenesis depends on the Majorana phases, so its understanding might provide the 
needed additional constraint. The baryon asymmetry can be expressed (Buchmuller et 
al 2002) in terms of the mass Mi of the lightest of the heavy Majorana neutrinos $i, 
the CP asymmetry e in <£>i decays, an effective neutrino mass rh, and the sum of the 
squares of the three light neutrino masses. Unfortunately, the relationship between e 
and the low-energy phases relevant for j3(5 is model dependent. Many models have been 
studied in the literature, each with its own relationship and conclusions. 

4. Calculating Nuclear Matrix Elements 

The observation of f3/3(0v) decay would immediately tell us that neutrinos are Majorana 
particles and give us an estimate of their overall mass scale but without accurate 
calculations of the nuclear matrix elements that determine the decay rate it will be 
difficult to reach quantitative conclusions about masses and hierarchies. Theorists have 
tried hard to develop many-body techniques that will allow such calculations. They have 
tried to calibrate their calculations to related observables: (3(3{2v) decay, ordinary (3 + 
and f3~ decay, Gamow- Teller strength distributions, odd-even mass differences, single- 
particle spectra. They have tried to exploit approximate isospin and SU (4) symmetries 
in the nuclear Hamiltonian, to extend well-known many-body methods in novel ways. 
In spite of all this effort, we know the matrix elements with limited accuracy. In this 
section we review the state of the nuclear-structure calculations and discuss ways to 
improve them. While increased accuracy is not easy to achieve, we will argue that it is 
also not impossible. 

Most recent attempts to calculate the nuclear matrix elements have been based on 
the neutron-proton Quasiparticle Random Phase Approximation (QRPA) or extensions 
to it. Of those that haven't, the most prominent are based on the shell model. While 
the two methods have much in common — their starting point is a Slater determinant 
of independent particles — the kinds of correlations they include are complementary. 
The QRPA treats a large fraction of the nucleons as "active" and allows these nucleons 
a large single-particle space to move in. But RPA correlations are of a specific and 
simple type best suited for collective motion. The shell model, by contrast, treats a 
small fraction of the nucleons in a limited single-particle space, but allows the nucleons 
there to correlate in arbitrary ways. That these very different approaches yield similar 
results indicates that both capture most of the important physics. 

4.1. QRPA 

The QRPA was developed by Halbleib and Sorenson (1967) and first applied to double- 
beta decay by Huffman (1970). Both it and early shell model calculations had problems 
reproducing measured /3/3(2v) rates until the realization by Vogel and Zirnbauer (1986) 
that in the QRPA the neutron-proton (np) particle-particle (i.e. pairing) interaction, 
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Figure 2. A consistency plot for the neutrino mass eigenvalues m\ and m 2 , for various 
hypothetical measurements. This set of curves indicates how measured values of S, 
{"mpp), <^ TO soi' an d ( m p) constrain the mass eigenvalues. The f3/3(0v) curve has been 
drawn for an incorrect value of the phase in the bottom panel to indicate the sensitivity 
of this technique for extracting the CP-violating phase. See text for further description 
of the parameters used to draw the curves. 
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Figure 3. A consistency plot for the neutrino mass eigenvalues mi, to 2 and 7713 for 
various hypothetical measurements. This set of surfaces indicates how measured values 
of E, (mpp), oscillations and (mp) constrain the mass eigenvalues. See text for further 
description of the parameters used to draw the surfaces. 

which has little effect on the collective Gamow- Teller resonance, suppresses /3(3{2v) 
rates considerably. Soon afterward, Engel et al (1988) and Tomoda and Faessler (1987) 
demonstrated a similar though smaller effect on fl/3(0u) decay. It was quickly realized, 
however, that the QRPA was not designed to handle realistic np pairing; the calculated 
half-lives were unnaturally sensitive to the strength of the pairing interaction. As a 
result, the rates of (3(3 decay, particularly (3{3(2i/) decay, were hard to predict precisely 
because a small change in a phenomenological parameter (the strength of np pairing) 
caused a large change in the lifetimes and eventually the breakdown (called a "collapse" ) 
of the entire method when the parameter exceeds some critical value. Most recent work 
in the QRPA has aimed at modifying the undesirable aspects of the method so that its 
sensitivity to np pairing becomes more realistic. 

One approach is "Second" QRPA (Raduta et al 1999) . The QRPA itself treats 
states in the intermediate nucleus as one-quasineutron + one quasiproton excitations 
of the initial and/or final ground states. The quasineutron-quasiproton pair is taken to 
obey boson statistics, inspiring the term "quasiboson approximation" for the method. 
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Second QRPA extends the structure of these states by adding another term in a boson 
expansion of the pair. 

A larger number of papers have been based on the Renormalized QRPA (RQRPA) 
(Hara 1964, Rowe 1968, and Suhonen 1995). The quasiboson approximation is 
equivalent to replacing commutators by their expectation values in an independent- 
quasiparticle approximation to the initial and final ground states. The RQRPA uses 
the QRPA ground states to evaluate the commutators. Because the commutators in 
turn help fix the ground states, the two must be evaluated self-consistently, usually 
via iteration. A variant of this approach is the "Full RQRPA" , in which the effects of 
isovector np pairing, artificially strengthened to account implicitly for isoscalar pairing, 
are included in the BCS calculation that defines the quasiparticles as well as in the 
subsequent QRPA calculation (Schwieger et al 1996; Simkovic et al 1997). [Isovector np 
pairing was first introduced in this way in the unrenormalized QRPA (Cheoun et al 1993, 
1995, Pantis et al 1996)]. Another extension is the Self-Consistent RQRPA (SCQRPA). 
Here, the BCS calculation is modified and iterated together with the RQRPA calculation 
until the RQRPA ground state has the same number of quasiparticles as the BCS-like 
state on which it is based. All these methods reduce the dependence of the (3(3(2v) 
matrix elements on the strength of the neutron-proton pairing interaction. 

Unfortunately, it's not clear which of the methods is best. One might presume that 
the Full RQRPA and the Self- Consistent RQRPA are better than the vanilla RQRPA, 
which in turn is better than the original unrenormalized QRPA. But the results of the 
RQRA and Full RQRPA appear to be quite sensitive to, e.g., the number of single- 
particle states in the model space (Simkovic et al 1997). They also violate an important 
sum rule for single (3 strength, and studies in solvable models suggest that the reduced 
dependence (at least of the RQRPA) on neutron-proton pairing may be spurious (Engel 
et al 1997) resulting from an artificial reduction of isoscalar pairing correlations. And it 
is not clear that the Full RQRPA's substitution of isovector pairing for isoscalar pairing 
is legitimate. 

The Self- Consistent QRPA has been applied to /3(3(0v) decay only once, and that 
application did not take into account the predictions of the same approach for (3(3{2v) 
decay. The approximations that go into the Second QRPA are of a different type, 
and while one study indicates less dependence on model space size for that method 
(Stoica and Klapdor-Kleingrothaus 2001) there are not many other reasons to prefer 
one approach or the other. 

Recently, the RQRPA was modified even further, so that the single (3 sum rule was 
restored. The resulting method, called the "Fully Renormalized QRPA" has yet to be 
applied to (3(3(01/) decay. Even more recently, Simkovic et al (2003) raised the issue of 
nuclear deformation, which has usually been ignored in QRPA-like treatments of nearly 
spherical nuclei^. The authors argued that differences in deformation between the initial 
and final nuclei can have large effects on the (3(3(2v) half-life. These ideas, too, have not 

% Psuedo-SU(3)-based truncations have been used to treat it in well-deformed nuclei (e.g., in Hirsch et 
al 1996). 
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yet been applied to (3(3(0v) decay. 

The profusion of RPA-based acronyms is both good and bad. The sheer number 
of methods applied gives us a kind of statistical sample of calculations, which we will 
use below to get an idea of the theoretical uncertainty in the matrix elements. But 
the sample may be biased by the omission of non-RPA correlations in all but a few 
calculations. Other approaches that include correlations more comprehensively should 
be pursued. 

4.2. Shell Model 

The obvious alternative to RPA, and the current method of choice for nuclear structure 
calculations in heavy nuclei where applicable, is the shell model. It has ability to 
represent the nuclear wave function to arbitrary accuracy, provided a large enough 
model space is used. This caveat is a huge one, however. Current computers allow very 
large bases (millions of states), but in heavy nuclei this is still not nearly enough. 
Techniques for constructing "effective" interactions and operators that give exact 
results in truncated model spaces exist but are hard to implement. Even in its crude 
form with relatively small model spaces and bare operators, however, the shell model 
offers advantages over the QRPA. Its complicated valence-shell correlations, which the 
QRPA omits (though it tries to compensate for them by renormalizating parameters) 
apparently affect the (3(3 matrix elements (Caurier et al 1996). 

The first modern shell-model calculations of (3(3 decay date from reference Haxton 
and Stephenson (1984) and references therein. Only a few truly large-scale shell model 
calculations have been performed. The heavy deformed (3(3 nuclei, 238 U, and 150 Nd, for 
example, require bases that are too large to expect real accuracy. Realistic work has 
thus been restricted to 48 Ca, 76 Ge, and 136 Xe, though less comprehensive calculations 
have been carried out in several other nuclei (Suhonen et al 1997). 

Large spaces challenge us not only through the problem of diagonalizing large 
matrices, but also by requiring us to construct a good effective interaction. The bare 
nucleon-nucleon interaction needs to be modified in truncated spaces (this is an issue in 
the QRPA as well, though a less serious one). Currently, effective interactions are built 
through a combination of perturbation theory, phenomenology, and painstaking fitting. 
The last of these, in particular, becomes increasingly difficult when millions of matrix 
elements are required. 

Related to the problem of the effective interaction is the renormalization of 
transition operators. Though the problem of the effective Gamow- Teller operator 
(Siiskonen et al 2001), which enters directly into (3(3{2v) decay, has drawn some 
attention, very little work has been done on the renormalization of the two-body 
operators that govern (3(3(0v) decay in the closure approximation. Shell model 
calculations won't be truly reliable until they address this issue, which is connected 
with deficiencies in the wave function caused by neglect of single-particle levels far from 
the Fermi surface. Engel and Vogel (2004) suggests that significant improvement on the 
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state of the art will be difficult but not impossible in the coming years. 
4-3. Constraining Matrix Elements with Other Observables 

The more observables a calculation can reproduce, the more trustworthy it becomes. 
And if the underlying model contains some free parameters, these observables can fix 
them. The renormalization of free parameters can make up for deficiencies in the model, 
reducing differences between, e.g., the QRPA and RQRPA once the parameters of both 
have been fit to relevant data. The more closely an observable resembles (3(3(Qv) decay, 
the more relevant it is. 

Gamow- Teller distributions, both in the (3~ and (3 + directions, enter indirectly 
into both kinds of (3(3 decay, and are measurable through (p,n) reactions. Aunola 
and Suhonen (1998) are particularly careful to reproduce those transitions as well as 
possible. Pion double charge exchange, in which a 7r + enters and a tt~ leaves, involves 
the transformation of two neutrons into two protons, like (3(3 decay, but the nuclear 
operators responsible aren't the same in the two cases. Perhaps the most relevant 
quantity for calibrating calculations of (3(3(0i/) decay is (3(3(2v) decay, which has now 
been measured in 10 different nuclei. 

Two recent papers have tried to use (3(3{2v) decay to fix the strength of np pairing 
in QRPA-based calculations. Stoica and Klapdor-Kleingrothaus (2001) used it only 
for the J n = 1 + channel relevant for (3(3{2v) decay, leaving the np pairing strength 
unrenormalized in other channels. By contrast, Rodin et al (2003) renormalized the 
strength in all channels by the same amount. The results of the two procedures were 
dramatically different: Stoica and Klapdor-Keingrothaus (2001) found that the /3/3(Oz/) 
matrix elements depended significantly on the theoretical approach (QRPA, RQRPA, 
FRQRPA, second QRPA) and, in some of the approaches, on the model space, while 
Rodin et al (2003) found almost no dependence on model-space size, on the form of 
the nucleon-nucleon interaction, or on whether the QRPA or RQRPA was used. The 
authors argued that fixing the np pairing strength to (3(3{2v) rates essentially eliminates 
uncertainty associated with variations in QRPA calculations of (3(3{Qv) rates, though 
they left open the question of how close to reality the calculated rates were. 

Given all these considerations, can we meaningfully estimate the uncertainty in the 
(3(3{^v) matrix elements? We address this question for a particular nucleus now. 

4-4- How Well Can We estimate Uncertainty? The Case of 7G Ge 

How accurate are existing calculations? One can answer provisionally by looking at the 
spread in the many predictions offered recently in the literature. We focus here on 76 Ge, 
which currently gives the strongest limit on (nipp) (or perhaps even a value for it: see 
discussions below), and figures prominently in several proposals for new experiments. 

Table EJ taken in part from Civitarese and Suhonen (2003), shows predictions of 
most of the calculations in the literature for the nuclear contribution to the decay rate 
(Doi et al 1985), C mm = (nipp) 2 T®Y 2 /ml = G 0u \M 0u \ 2 ml, and for the effective 
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neutrino mass (mpp) that would be deduced by measuring a lifetime Ty 2 = 4.0 x 10 27 
years. Besides differing in method (by which they are grouped) the calculations 
use different model spaces, fit different observables, and adjust different parameters. 
The resulting C mm 's vary by nearly two orders of magnitude, leading to an order of 
magnitude variation in the extracted neutrino mass. But some of the "outliers", labeled 
with asterisks in the table, can be eliminated on the grounds that the corresponding 
calculations omit or misrepresent clearly important effects. The VAMPIR model used 
in the early calculation of Tomoda et al (1986) contains no np pairing correlations and 
so gives too large a decay rate. The shell-model truncation by Haxton and Stephenson 
(1984) was done in a way that minimized np pairing, and the upper limit in C mm from 
Engel et al (1989) was considered probably too large by the authors (though in their 
calculation they set qa — 1, which here would move the upper limit into the middle of 
the range). At the other end of the spectrum, the very small self-consistent RQRPA 
decay rates of reference Bobyk et al (2001) were obtained with a value of the np pairing 
strength that was not consistent with the measured (3(3{2v) rate; when the strength is 
adjusted to reproduce (3(3(2v) decay, the results for the PP(0u) rate are close to those 
of the plain QRPA in the same reference. Without any further culling, the remaining 
C's vary by about 1 order of magnitude, and the extracted (m^'s vary by a factor of 
about three, from 0.022 to 0.068 eV for the lifetime we have chosen. 

Even if some of the other calculations are objectionable, it is difficult to reduce the 
spread much below this factor of three without some real work. Aunola and Suhonen 
(1998), who among all the entries do the most extensive job of adjusting parameters to 
reproduce spectroscopic data and single-/? transition strengths, obtain a low neutrino 
mass: 0.022 eV from our hypothetical lifetime. Caurier et al (1996), who do the most 
complete shell-model calculation, obtain 0.059 eV, close to the maximum mass from 
calculations we haven't removed. We cannot discount these two careful calculations 
without examining them more closely and so cannot further restrict the range. This is 
not to say that these calculations can't be questioned. It may not be appropriate, for 
example, to fit single-particle energies to /3-decay rates as Aunola and Suhnen do. But 
for now, until someone investigates the situation further, the range of reasonable results 
should include these two calculations. 

As mentioned above, Rodin et al (2003) argue that the variation in QRPA and 
RQRPA rates can be nearly eliminated by renormalizing a few parameters to reproduce 
pairing gaps and (3{3{2v) decay rates. Indeed, the table shows that the QRPA and 
RQRPA numbers from that reference are close. But it has yet to be shown conclusively 
that (3/3 (2v) decay is more important than single-/? decay or other observables, and 
Engel and Vogel (2004) demonstrate that in a two- level SO (5)-based model (Dussel et 
al 1970), at least, the procedure cannot fully eliminate model-space dependence when 
the QRPA is used. 
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Table 2. Predictions for 76 Ge. The quantity C mm , defined in the text, is a measure of 
the nuclear contribution to the f3f3(0v) decay rate. The values of {mpp) are what those 
that would be extracted by a calculation if the lifetime were 4.0 x 1— 27 years. The 
asterisks indicate approaches that omit important physics or have been superceded. 
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4-5. Reducing the Uncertainty 

What can be done to improve the situation? In the near term, improvements can be 
made in both QRPA-based and shell- model calculations. First, existing calculations 
should be reexamined to check for consistency. One important issue is the proper 
value of the axial-vector coupling constant qa, which is often set to 1 (versus its real 
value of 1.26) in calculations of (3 decay and f3(3(2u) decay to account for the observed 
quenching of low-energy Gamow- Teller strength. What value should one use for PP(0u) 
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decay, which goes through intermediate states of all multipolarity, not just 1 + ? Some 
authors use Qa — 1, some qa = 1.26, and some Qa = 1 for the 1 + multipole and 1.26 
for the others. (Often authors don't reveal their prescriptions.) The second of these 
choices appears inconsistent with the treatment of f3(3(2v) decay. Since the square of 
Qa enters the matrix element, this issue is not totally trivial. The striking results of 
Rodin et al suggest that an inconsistent treatment is responsible for some of the spread 
in Table 121 More and better charge-exchange experiments would help teach us whether 
higher-multipole strength is also quenched. 

Next, the various versions of the QRPA should be tested against exact solutions in 
a solvable model that is as realistic as possible. The most realistic used so far are the 
S'0(5)-based model first presented by Dussel et al (1970) and used to study the QRPA 
and RQRPA for Fermi (3j3{2v) decay by Hirsch et al (1997), a two-level version of that 
model used by Engel and Vogel (2004) for the QRPA in Fermi /3/3(2z/) and /3/3(0z/) decay, 
and an S'0(8)-based model (Evans et al 1981) used to test the QRPA and RQRPA 
for both Fermi and Gamow- Teller f3f3(2v) decay by Engel et al (1997). It should be 
possible to extend the 50(8) model to several sets of levels and develop techniques for 
evaluating j3f3(0v) matrix elements in it. All these models, however, leave out spin-orbit 
splitting, which weakens the collectivity of np pairing. A generalization of pseudospin- 
based models like those of Ginocchio (1980), used in the Fermion Dynamical Symmetry 
Model (Wu et al 1994), might be useful provided techniques like those of Hecht (1993) 
to evaluate expectation values of operators outside the algebra can be extended. Such 
calculations should help us understand the virtues and deficiencies of QRPA extensions. 

Along the same lines, we will need to understand the extent to which such methods 
can reproduce other observables, and their sensitivity to remaining uncertainties in 
their parameters. A very careful study of the first issue was made by Aunola and 
Suhonen (1998) and the second has been explored in several papers, most thoroughly 
by Stoica and Klapdor-Kleingrothaus (2001). These efforts must be extended. The 
work is painstaking, and perhaps not as much fun as concocting still more variations of 
the QRPA, but it is crucial if we are to reduce theoretical uncertainty. Self-consistent 
Skyrme HFB+QRPA, applied to single-/? decay by Engel et al (1999) and Gamow- 
Teller resonances by Bender et al (2002), may be helpful here; it offers a more general 
framework, in principal anyway, for addressing the variability of calculated matrix 
elements. Solvable models can be useful here too, because they can sometimes supply 
synthetic data to which parameters can be adjusted (as in Engel and Vogel (2004)). 

The best existing shell-model calculation produces smaller matrix elements than 
most QRPA calculations. Computer speed and memory is now at the point where 
the state of the shell-model art can be improved. The calculation of the f}(3 decay 
of 76 Ge by Caurier et al (1996) used the /5/2P3/2Pi/25 , 9/2 model space, allowing up to 
8 particles (out of a maximum of 14) into the g§/2 level. Nowadays, with the help 
of the factorization method (Papenbrock and Dean 2003; Papenbrock et al 2003a), an 
accurate approximation to full shell-model calculations, we should be able to fully occupy 
the gg/2 level, and perhaps include the g 7 / 2 and / 7 / 2 levels (though those complicate 
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things by introducing spurious center-of-mass motion). In addition, one can try through 
diagrammatic perturbation theory to construct effective (3/3(Qv) operators for the model 
space that are consistent with the effective interaction. Though perturbation theory has 
convergence problems, the procedure should at least give us an idea of the uncertainty 
in the final answers, perhaps also indicating whether result obtained from the "bare" 
operators is too large or too small. Research into effective operators has been revived 
in recent years (Haxton and Song 2000) and we can hope to improve on diagrammatic 
perturbation theory. One minor source of uncertainty connected with renormalization 
(which also affects the QRPA) is short-range two-nucleon correlations, currently treated 
phenomenologically, following Miller and Spencer (1976). 

Some ab initio calculations of nuclear properties are now possible (Carlson 1998). 
Although researchers have obtained accurate results only in nuclei with A < 12, we have 
reason to hope that "nearly exact" variational or coupled-cluster methods (Kowalski et 
al 2003) will be applied to medium-heavy nuclei such as 76 Ge in the next 5 years. 
Accurate calculations in heavier nuclei are probably further off. 

In short, much can be done and we would be well served by coordinated attacks on 
these problems. There are relatively few theorists working in (3(3 decay, and their efforts 
have been fragmented. More collaborations, postdoctoral and Ph.D projects, meetings, 
etc., would make progress faster. There is reason to be hopeful that the uncertainty will be 
reduced. The shell-model matrix element may be too small because it does not include 
any particles outside the /p-shell. These particles, as shown by QRPA calculations, 
make the matrix element larger. We suspect that the results of a better shell-model 
calculation will be closer than the best current one to the QRPA results and that, as 
noted above, the spread in those results can be reduced. Finally, other nuclei may be 
more amenable to a good shell-model calculation than Ge. 136 Xe has 82 neutrons (a 
magic number) making it a particularly good candidate. 

5. Other Possible Mechanisms for Double-Beta Decay 

Although the occurrence of (3(3(0u) decay implies the existence of massive Majorana 
neutrinos (Schechter and Valle 1982), their exchange need not be the dominant 
contribution to the decay rate. Almost any physics that violates lepton number can 
cause (3(3{Qv) decay. A heavy Majorana neutrino can be exchanged, or supersymmetric 
particles, or a leptoquark. Right-handed weak currents, either leptonic or hadronic, 
can cause the absorption of an emitted virtual neutrino without the helicity flip that 
depends on the neutrino mass. These possibilities have been reviewed by Faessler and 
Simkovic (1998) and Suhonen and Civitarese (1998). But now we know that there are 
light neutrinos and that next-generation (3(3{^v) experiments may well allow us to learn 
something new about them. The other possibilities are more speculative and instead of 
analyzing them in detail we confine ourselves to the question of whether it is possible, 
should P(3(0u) decay be observed, to determine which mechanism is responsible. 

Light-neutrino exchange with right-handed currents is unique in this regard because 
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the contraction of left and right-handed currents gives rise to a term g p 7 p in the 
numerator of the neutrino propagator (see Equation [SJ) that cancels in the contraction 
of two left-handed or two right-handed currents. The extra q allows the electron to be 
emitted in a p-wave as well as an s-wave and introduces a contribution to the amplitude 
from nucleon recoil. In addition to increasing the phase space, these effects lead to 
a different single-electron energy distribution and opening-angle dependence than in 
/3/3(Ozv) decay driven by the neutrino mass. An experiment sufficiently sensitive to 
the energies and paths of individual electrons could therefore determine whether right- 
handed currents were driving the decay. 

The effects of right-handed currents can be shown to be negligible unless there 
exist some neutrinos with masses larger than a few MeV. Heavy Majorana neutrinos, 
however, are generated naturally in a variety of models via the see-saw mechanism. The 
heavy neutrinos complicate matters because they can mediate /3/3(0z/) decay themselves. 
It seems possible within, e.g., left-right symmetric models, for the exchange of these 
neutrinos via right-handed currents to compete with or even dominate the exchange of 
light-neutrinos. 

The exchange of heavy particles involves short-range propagators that give rise to 
decay rates of the same form as in the mass-driven mode: simple two-s-wave-electron 
phase space multiplied by the square of an amplitude§. The angular distributions and 
single-electrons spectra will therefore be the same in all these processes. The only way to 
distinguish one from another is to take advantage of the different nuclear matrix elements 
that enter the amplitudes (leading to different total decay rates). Unknown parameters 
such as the effective neutrino mass or the trilinear i?-parity-violating supersymmetric 
coupling (violation of R parity naturally accompanies Majorana neutrino-mass terms) 
also enter the rates, so several transitions would have to be measured. Simkovic and 
Faessler (2002) argue that this is best accomplished by measuring transitions to several 
final states in the same nucleus, but if the matrix elements can be calculated accurately 
enough one could also measure the rates to the ground states of several different nuclei. 

The problems in determining the source of /3/5(Oz/) decay are mitigated by 
constraints from other experiments on many extra-Standard models. Some of these 
constraints will be much stronger once the Large Hadron Collider comes on line. If no 
signs of supersymmetry (for example) appear there, then supersymmetry probably does 
not exist. If supersymmetric particles are found, but neutralinos make it out of the 
detector without decaying, then .R-parity, which prevents the decay of supersymmetric 
partners into familiar particles, will not be strongly violated. The trilinear i?-parity- 
violating coupling could then be ruled out as the source of double-beta decay. The 
presence of the LHC will make complementary experiments on (3(3(0u) decay still more 
attractive. 

Some theories posit the emission Goldstone bosons called "Majorons" , together 
with the emission of electrons, in /5/3(Oz/) decay. If Majorons were emitted in a detector 

§ Two-nucleon correlations do not suppress the effects of heavy particles, which can be transmitted 
between nucleons by pions (Prezeau et al 2003, Faessler et al 1997). 
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the total energy carried by the newly created electrons would vary from event to event. 
A long lifetime for this decay would make it difficult to detect above a /5/5(2z/)-decay 
background. Double-beta decay has also been discussed as as test of special relativity 
and the equivalence principle (Klapdor-Kleingrothaus et al 1999). Finally, very recent 
attempts to unify the dark sector and neutrino physics (Kaplan et al 2004) posit a scalar 
field with variations on the scale of millimeters that couples to neutrinos. If such a field 
existed, the rate of (3(3(0u) decay might depend on the density of matter in which the 
process occurred. 

6. Experimental Situation 

If an experiment observes /3/3(0i/) it will have profound physics implications. Such an 
extraordinary claim will require extraordinary evidence. The recent claim (Klapdor- 
Kleingrothaus et al 2004) for an observation of (3f3(0u) has been controversial (See 
discussion below). Also previous "false peaks" in (3(3 spectra have appeared near a 
(3(3(0v) endpoint energy (see discussion in Moe and Vogel 1994, page 273). One must 
ask the question: What evidence is required to convincingly demonstrate that (3(3(Qv) 
has been observed? Low-statistical-significance peaks (« 2a) have faded with additional 
data, so one must require strong statistical significance (perhaps 5a). (See Fig. HJ) This 
will require a large signal-to-noise ratio that will most likely be accomplished by an ultra- 
low-background experiment whose source is its detector. Such experiments are usually 
calorimetric and provide little information beyond just the energy measurement. 

How does an experiment demonstrate that an observed peak is actually due to (3(3 
decay and not some unknown radioactivity? Additional information beyond just an 
energy measurement may be required. For example, although there is some uncertainty 
associated with the matrix elements, it is not so large that a comparison of measured 
rates in two different isotopes could not be used to demonstrate consistency with the 
Majorana-neutrino hypothesis. Alternatively, experiments that provide an additional 
handle on the signal, for example by measuring a variety of kinematical variables, 
demonstrating that 2 electrons are present in the final state, observing the 7 rays 
associated with an excited state, or identifying the daughter nucleus, may lend further 
credibility to a claim. Experiments that provide this extra handle may require a 
significantly more complicated apparatus and therefore face additional challenges. 

The exciting aspect of (3(3 research today is that many proposed experiments 
intend to reach a Majorana mass sensitivity of \J o~rn 2 Atm . Several different isotopes and 
experimental techniques are being pursued actively and many of the programs look 
viable. In this section we describe the current situation in experimental (3(3 decay . 

6.1. Results to Date 

Table [HI lists the recent (3(3{0i/) results. The best limits to date come from the enriched 
Ge experiments. The two experiments had comparable results although the Heidelberg- 
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Figure 4. This contour plot shows the half life, in units of 10 25 y, for a peak of 5<r 
significance for a given exposure and background. The KDHK point is shown. 

Moscow result was marginally better. The TfY 2 limits near 2 x 10 25 y results in a {mpp) 
limit near 300 meV, with an uncertainty of about a factor of 3 because of the uncertainty 
in |Moj/|. One recent paper (Zdesenko et al 2002) performed a joint analysis of the two 
experiments and found Ty 2 > 2.5 x 10 25 y. 

Most of the results listed in Table El are at least a few years old. The obvious 
exceptions to this are the Te and Cd results. CUORICINO continues to collect data. 

6.1.1. A Claim for the Observation of /3/3(0u) In early 2002, a claim for the observation 
of PP(0i>) was published (Klapdor-Kleingrothaus et al 2002a). The paper made a 
poor case for the claim and drew strong criticism (Aalseth et al 2002a, Feruglio et 
al 2002, Zdesenko et al 2002). The initial response to the criticism was emotional 
(Klapdor-Kleingrothaus 2002b). In addition, one of the original co-authors wrote a 
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Table 3. A summary of the recent f3P(0v) results. The (mpp) limits are those deduced 
by the authors. All limits are at 90% confidence level unless otherwise indicated. The 
columns providing the exposure and background are based on arithmetic done by the 
authors of this paper, who take responsibility for any errors in interpreting data from 
the original sources. 



Isotope 


Exposure 


Background 


Half-Life 






(kmole-y) 


(counts) 


Limit (y) 


(meV) 


48 Ca 


5 x KT 5 





> 1.4 x 10 22 


< 7200 - 44700 c 


76 Ge 


0.467 


21 


> 1.9 x 10 25 


< 350 fe 


76 Ge 


0.117 


3.5 


> 1.6 x 10 25 


< 330 - 1350 c 


76 Ge 


0.943 


61 


= 1.2 x 10 25 


= 440 d 


82 Se 


7 x 10~ 5 





> 2.7 x 10 22 (68%) 


< 5000 e 


100 Mo 


5 x 10~ 4 


4 


> 5.5 x 10 22 


< 2100 / 


116 Cd 


1 x lO" 3 


14 


> 1.7 x 10 23 


< noo 9 


128 Te 


Geochem. 


NA 


> 7.7 x 10 24 


< 1100 - 1500 h 


i30 Te 


0.025 


5 


> 5.5 x 10 23 


< 370 - 1900* 


136 Xe 


7 x 10~ 3 


16 


> 4.4 x 10 23 


< 1800 - 5200 J ' 


i50 Nd 


6 x 10~ 5 





> 1.2 x 10 21 


< 3000 fc 



a Ogawa et al 2004; b Klapdor-Kleingrothaus et al 2001; c Aalseth et al 2002; d Klapdor- 
Kleingrothaus et al 2004; e Elliott et al 1992; / Ejiri et al 2001; 9 Danevich et al 2003; 
h Bernatowicz et al 1993; l Arnaboldi et al 2004; 3 Luescher et al 1998; fc De Silva et al 
1997 

separate reply (Harney 2001) that mostly defended the claim yet acknowledged some 
significant difficulty with the analysis. This author's name doesn't appear on later 
papers. More recently, however, supporting evidence for the claim has been presented 
and we recommend the reader study Klapdor-Kleingrothaus et al (2002c) for a good 
discussion of the initial evidence and Klapdor-Kleingrothaus et al (2004) for the most 
recent data analysis. Importantly, this later paper includes additional data and therefore 
an increase in the statistics of the claim. In this subsection we summarize the current 
situation. (We use the shorthand KDHK to refer to the collection of papers supporting 
the claim.) 

FigureElshows the spectrum corresponding to 71.7 kg-y of data from the Heidelberg- 
Moscow experiment between 2000 and 2060 keV (Klapdor-Kleingrothaus et al 2004). 
This spectrum is shown here to assist the casual reader in understanding the issues. 
However, the critical reader is encouraged to read the papers listed in the references as 
the authors analyze several variations of this data using different techniques. The fit 
about the expected /3/3(0z/) peak energy yields 28.75 ± 6.86 counts assigned to (3f3{Qp). 
The paper claims a significance of approxmately 4a for the peak, where the precise 
significance value depends on the details of the analysis. The corresponding best-fit 
lifetime, Ty 2 = 1.19 x 10 25 years (Klapdor-Kleingrothaus et al 2004), leads to a (m^) 
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Figure 5. The spectrum from the Heidelberg- Moscow experiment upon which the 
claim for (3(3(Qv) is based. The data in the two panels are identical. The lower panel 
has a Gaussian curve to indicate the strength of the claimed [}[3{§v) peak. 



of 440 meV with the matrix element calculation of Staudt et al (1990) chosen by the 
authors. 

In the region between 2000 and 2100 keV, the KDHK analysis of 2002 found a total 
of 7 peaks. Four of these were attributed to 214 Bi (2011, 2017, 2022, 2053 keV), one 
was attributed to (3(3(0v) decay (2039 keV), and two were unidentified (2066 and 2075 
keV). The KDHK analysis of 2004 does not discuss the spectrum above 2060 in detail. 
An additional possible feature may also be present near 2030 keV. A study (Klapdor- 
Kleingrothaus et al 2003a) comparing simulation to calibration with 214 Bi demonstrates 
that if the location of the Bi is known, the spectrum can be calculated. Furthermore, 
the relative strengths of the strong Bi lines at 609, 1764 and 2204 keV can be used to 
determine the location of the activity. Because the results of summing effects depend 
on the proximity of the activity, its location is critical for the simulation of the weak 
peaks near the (3(3(0v) endpoint. The study also shows that the spectrum can't be 
naively estimated, as was done in Aalseth et al (2002a). In fact, Table VII in Klapdor- 
Kleingrothaus (2002c) finds, even with a careful simulation, that the expected strengths 
of the 214 Bi peaks in the 2000-2100 keV region are not predicted well by scaling to 



CONTENTS 



28 



the strong peaks. That is, the measured intensities of the weak peaks are difficult to 
simulate without knowing the exact location of the activity. Furthermore, the deduced 
strengths of the weak lines are more intense than expected by scaling from the strong 
peaks, even though the activity location is chosen to best describe the relative intensities 
of the strong peaks. 

Double-beta decay produces two electrons that have a short range in solid Ge. 
Therefore, the energy deposit is inherently localized. Background process, such as the 
7 rays from 214 Bi, tend to produce multiple energy deposits. The pulse waveform can 
be analyzed to distinguish single site events (SSE) from multiple site events. Such an 
analysis by KDHK (Klapdor-Kleingrothaus et al 2003c, 2004) tends to indicate that the 
Bi lines and the unidentified lines behave as multiple site events, whereas the (3(3(0v) 
candidate events behave as SSE. Note, however, that the statistics are still poor for the 
experimental lines and this conclusion has a large uncertainty. Nonetheless, this feature 
of the data is very intriguing and clearly a strength of the KDHK analysis. 

An analysis by Zdesenko et al (2003) points out the strong dependence of the result 
on the choice of the window width of the earlier 2002 analysis. The KDHK analysis 
argues that a small window is required because of the neighboring background lines. 
Even so, their Monte Carlo analysis shows that the result becomes less stable for small 
windows (see Fig. 9 in Klapdor-Kleingrothaus et al 2002c). Zdesenko et al (2003) 
also remind us that the significance of a signal is overestimated when the regions used 
to estimate the background are comparable to the region used to determine the signal 
(Narsky 2000). The report of Klapdor-Kleingrothaus et al (2004) fits a wide region 
containing several peaks simultaneously after using a Bayesian procedure to identify the 
location of the peaks. 

The claim for f3f3(0v) decay was made by a fraction of the Heidelberg-Moscow 
collaboration. A separate group of the original collaboration presented their analysis 
of the data at the IV International Conference on Non-Accelerator New Physics 
(Bakalyarov et al 2003). They indicate that the data can be separated into two distinct 
sets with different experimental conditions. One set includes events that are described as 
"underthreshold pulses" and one set that does not. Analysis of the two sets produce very 
different conclusions about the presence of the claimed peak. They conclude that the 
evidence is an experimental artifact and not a result of /9/9(0z/)decay. KDHK responds 
(Klapdor-Kleingrothaus et al 2004) that these corrupt data were not included in their 
analysis. 

Traditionally, (3(3 experiments have ignored systematic uncertainties in their 
analysis. Only recently with the start-up of high-statistics (3(3(2v) results has this 
situation begun to change. Historically, /3(3(0v) results have always been quoted as 
upper limits based on low count rates. As a result, systematic uncertainties tended to be 
negligible in the final quoted values. With a claim of a positive result, however, the stakes 
are dramatically raised. It is clear that it is difficult to produce a convincing result when 
the signal counts are comparable to expected statistical fluctuations in the background. 
The further presence of nearby unidentified peaks makes the case even harder to prove. 
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Although KDHK does discuss some systematic uncertainties qualitatively and indicates 
they are small (in the position of the (3(3(0v) peak, and the expected peak width, for 
example), there is no consideration of an uncertainty associated with the background 
model. 



The next round of proposed (3(3(0u) experiments are designed to reach J 5ml, 



and therefore will quickly confirm or repudiate this claim. This is fortunate since the 
feature near 2039 keV in the KDHK claim will likely require an experimental test. These 
experiments should provide a detailed listing of all identified systematic uncertainties 
and a quantified estimate of their size. Furthermore, because the stakes are very high 
and there will be many people who are biased, either for or against the KDHK claim, 
blind analyses should also become part of the experimental design. 

6.2. Future Possibilities for (3 /3(0^) Experiments 

The recent review by Elliott and Vogel (2002) describes the basics of experimental 
(3(3(0i>) decay in some detail. Therefore, we refer the reader to that article and only 
summarize the status of the various projects. Table |H lists the proposals. 

6.2.1. CANDLES The CANDLES collaboration has recently published the best limit 
on PP(0p) decay of 1.4 x 10 22 y in 48 Ca (Ogawa et al 2004). Using the ELEGANTS VI 
detector, this experiment consisted of 6.66 kg of CaF 2 (Eu) crystals surrounded by Csl 
crystals, a layer of Cd, a layer of Pb, a layer of Cu, and a layer of LiH-loaded paraffin, all 
enclosed within an air-tight box. This box was then surrounded by boron-loaded water 
tanks and situated underground at the Oto Cosmo Observatory. This measurement 
successfully demonstrated the use of these crystals for (3(3 studies. 

An improved version of this crystal technology, the CANDLES-III detector 
(Kishimoto et al 2004), is presently being constructed with 200 kg of CaF 2 crystals. 
These crystals have better light transmission than the CaF 2 (Eu) crystals. This design 
uses sixty 10-cm 3 CaF 2 crystals, which are immersed in liquid scintillator. The 
collaboration has also proposed a 3.2-t experiment that hopes to reach 100 meV for 

6.2.2. COBRA The COBRA experiment (Zuber 2001) uses CdZnTe or CdTe 
semiconductor crystals. These crystals have many of the advantages of Ge detectors 
but, in addition, operate at room temperature. Because the crystals contain Cd and 
Te, there are 7 (3(3 and (3 + (3 + isotopes contained. The final proposed configuration is 
for 64000 1-cm 3 crystals for a total mass of 370 kg. The collaboration has already 
obtained 30-keV resolution at 2.6 MeV with these detectors and has published initial 
/3/3-decay studies (Kiel et al 2003). Background studies are the current focus of the 
efforts. Although it is tempting to focus on the naturally isotopic abundant 130 Te for 
(3{3(0v) decay, the presence of the higher Q-value 116 Cd creates a serious background 
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Table 4. A summary of the (3(3{Qv) proposals. Background estimates were not 
available for all projects. The quantity of isotope includes the estimated efficiency 
for fi/3(0v). 



Collaboration 


Isotope 


Anticipated 


Detector 




(kmol) 


Background 
(counts/y) 


Description 


CAMEO" 


116 Cd (2) 


few/year 


CdW04 crystals in liq. scint. 


CANDLES 6 


48 Ca (0.04) 




CaF 2 crystals in liq. scint. 


COBRA c 






CdTe semiconductors 


CUORE rf 


130 Te (1.4) 


~ 60/y 


Te02 bolometers 


DCBA e 


82 Se (2) 


~ 40/y 


Nd foils and tracking chambers 


EXO / 


136 Xe (4.2) 


<i/y 


Xe TPC, 


GEM 9 


76 Ge (11) 


~ 0.8/y 


Ge detectors in LN 


GENIUS' 1 


76 Ge (8.8) 


~ 0.6/y 


Ge detectors in LN 


GSO* 


160 Gd (1.7) 




Gd2Si05 crystals in liq. scint. 


Majorana J 


76 Ge (3.5) 


~ 1/y 


Segmented Ge detectors 


MOON fc 


100 Mo (2.5) 


~8/y 


Mo foils and plastic scint. 


MPI bare Ge' 


76 Ge (8.8) 




Ge detectors in LN 


nano- crystals™ 


100 kmol 




suspended nanoparticles 


super-NEMO™ 


82 Se (0.6) 


~ 1/y 


foils with tracking 


Xe° 


136 Xe (6.3) 


~ 118/y 


Xe dissolved in liq. scint. 


XMASS P 


136 Xe (6.1) 




liquid Xe 



a Bellini et al 2001; fe Kishimoto et al 2004; c Zuber 2001; d Arnaboldi et al 2004a; e Ishihara 
2000; -fDanilov et al 2000; ^Zdesenko et al 2001; ^Klapdor-Kleingrothaus et al 2001; 
^Danevich et al 2001, Wang et al 2000; ^'Gaitskell et al 2003; fc Ejiri et al 2000; 'Abt 
et al 2004; "McDonald 2004; "Sarazin et al 2000; °Caccianiga and Giammarchi 2001; 
p Moriyama et al 2001 

from its (3(3(2v) decay. Detectors enriched in 116 Cd are probably required to reach 45 
meV. 

6.2.3. CUORE The CUORICINO experiment uses 41 kg of Te0 2 crystals operated at 
10 mK as bolometers. During the initial cool down, some of the cabling failed and hence 
not all crystals were active. As a result the initial run had contained about 10 kg of 130 Te 
(Arnaboldi et al 2004). An initial exposure of 5.46 kg-y, with an energy resolution of 
9.2 keV FWHM resulted in a limit Tfy 2 > 7.2 x 10 23 y at 90% confidence level (Norman 
2004). The background in the region of interest for this run was 0.22 ±0.04 counts/ (keV 
kg y). The experiment has been suspended in order to fix the cabling. Afterward, a 
3-year run with the full mass will have a sensitivity of 10 25 years. 

The CUORICINO project is a prototype for the CUORE proposal. CUORE would 
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contain 760 kg of TeCV With the anticipated improvement in background to better 
than 0.01 counts/ (keV kg y), the half life sensitivity is ~ 7 x 10 26 y or a few 10's of meV 
for (rapp) (Arnaboldi et al 2004a). 

6.2.4. DCBA The Drift Chamber Beta-ray Analyzer (DCBA) (Ishihara 2000) is a 
tracking chamber within a 1.6-kG magnetic field that can examine any (3(3 source that 
can be formed into a thin foil. On both sides of the foil are tracking regions filled with 
1-atm He gas. The future program calls for 25 m 2 of source foils contained within each 
of 40 modules, for a total of approximately 600 kg of source material. The half-life 
sensitivity is estimated to be a few xlO 26 y for 82 Se, 100 Mo, or 150 Nd, assuming an 
enrichment of 90%. 

6.2.5. EXO The EXO project proposes to use 1-10 t of about 80% enriched liquid Xe 
as a time projection chamber (Danilov et al 2000). Development of a high-pressure gas 
TPC is being pursued in parallel. In addition to measuring the energy deposit of the 
electrons, the collaboration is developing a technique for extracting the daughter Ba ion 
from the Xe and detecting it offline. Observing the daughter in real time with the (3(3 
decay is a powerful technique for reducing background. With a 1-ton experiment, they 
anticipate sensitivity to a lifetime of 8 x 10 26 y. 

The collaboration has had some good progress on the research and development 
required to demonstrate that this technically challenging project is feasible. They 
have determined the energy resolution by using both ionization and scintillation 
measurements in liquid Xe. The resolution result a = 3% stated in Conti et al (2003) 
was measured at 570 keV. Assuming a statistical dependence on energy this means 
about 1.5% resolution at the (3(3(Qv) energy of 2480 keV. They have also built an atom 
trapping system and have observed lone Ba ions in an optical trap. Furthermore, they 
have begun experiments to demonstrate that the ions are trapped and observable in 
an appreciable Xe gas background (Piepke 2004). Finally, using a 222 Ra source they 
are testing the Ba extraction technology. Ra and Ba have similar chemistry, but the 
radioactive decay of Ra makes it a convenient test material. 

The EXO team is currently preparing a 200-kg enriched- Xe experiment to operate 
at the Waste Isolation Pilot Plant (WIPP). This prototype will not initially include Ba 
extraction. 

6.2.6. MOON The MOON project (Ejiri et al 2000) proposes to use 1 1 of Mo enriched 
to 85% in 100 Mo. The beauty of 100 Mo is that it not only is a good /3(3(0u) isotope, 
but also has a large charge-current cross section for low-energy solar neutrinos. Thus 
the MOON detector is being designed to perform both experiments. MOON measures 
individual (3 rays from (3(3 decay, which helps identify events arising from the light- 
neutrino-mass (3(3 mechanism and improve the background rejection. 

The reference design for MOON is a collection of modules of interleaved plate and 
fiber scintillators sandwiching Mo foils. Each foil is about 20 mg/cm 2 . Good position 
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resolution is required to exploit the timing of the radioactive product produced in the 
solar-neutrino interaction. The position is determined by the fiber scintillators, whereas 
the scintillator plate provides the energy resolution (a ~ 2.2% at 3 MeV). Other detector 
options are under consideration but a 1-kg prototype of the reference design is currently 
being built. 

6.2. 7. Majorana The Majorana Collaboration proposes to field 500 kg of 86% enriched 
Ge detectors (Gaitskell et al 2003). By using segmented crystals and pulse-shape 
analysis, multiple-site events can be identified and removed from the data stream. 
Internal backgrounds from cosmogenic radioactivities will be greatly reduced by these 
cuts and external 7-ray backgrounds will also be preferentially eliminated. Remaining 
will be single-site events like that due to (3(3. The sensitivity is anticipated to be 4 x 10 27 

y- 

Several research and development activities are currently proceeding. The 
collaboration is building a multiple-Ge detector array, referred to as MEGA, that will 
operate underground at the Waste Isolation Pilot Plant (WIPP) near Carlsbad, NM 
USA. This experiment will investigate the cryogenic cooling of many detectors sharing a 
cryostat in addition to permitting studies of detector-to-detector coincidence techniques 
for background and signal identification. A number of segmented crystals are also being 
studied to understand the impact of segmentation on background and signal. This SEGA 
program consists of one 12-segment enriched detector and a number of commercially 
available segmented detectors. Presently, commercially available segmented detectors 
are fabricated from n-type crystals. Such crystals are much more prone to surface 
damage and thus more difficult to handle when packaging inside their low-background 
cryostats. Hence the collaboration is also experimenting with segmenting p-type 
detectors. 

6.2.8. Bare Ge Crystals The GENIUS collaboration (Klapdor-Kleingrothaus et al 
2001) proposed to install 1 t of enriched bare Ge crystals in liquid nitrogen. 
By eliminating much of the support material surrounding the crystals in previous 
experiments, this design is intended to reduce backgrounds of external origin. Note 
how this differs from the background-reduction philosophy associated with pulse-shape 
analysis coupled with crystal segmentation. The primary advocates for this project 
indicate (Klapdor-Kleingrothaus et al 2004) that its motivation has been questioned 
by their own claim of evidence for /3/3(Oz/)decay. Even so, the GENIUS test facility 
(Klapdor-Kleingrothaus et al 2003) is being operated to demonstrate the effectiveness 
of operating crystals naked in liquid cryogen. 

Another group at the Max Plank Institute in Heidelberg, however, is proposing to 
pursue a similar idea. They have recently submitted a Letter of Intent (Abt 2004) to the 
Gran Sasso Laboratory. They propose to collect the enriched Ge crystals from both the 
Heidelberg-Moscow and IGEX expeiments and operate them in either liquid nitrogen or 
liquid argon. As a second phase of the proposal, they plan to purchase an additional 20- 
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kg of enriched Ge detectors (most likely segmented) and operate with a total of 35 kg for 
about 3 years. Finally, they eventually plan to propose a large ton-scale experiment. It 
should be noted that this collaboration and the Majorana collaboration are cooperating 
on technical developments and if a future ton-scale experiment using 76 Ge proceeds 
these two groups will most likely merge and optimally combine the complementary 
technologies of bare-crystal operation and PSA-segmentation. 

6.3. Nanocrystals 

Some elements may be suitable for loading liquid scintillator with metallic-oxide 
nanoparticles. Since Rayleigh scattering varies as the sixth power of the particle radius, 
it can be made relatively small for nanoparticles of radii below 5 nanometers. Particles 
of this size have been developed and commercial suppliers of Zr0 2 , Nd 2 3 etc. are 
available. Absorption of the materials must also be taken into account, but some of the 
metal oxides such as Zr0 2 and Te0 2 are quite transparent in the optical region because of 
the substantial band gaps in these insulators. Some members of the SNO collaboration 
(McDonald 2004) have been studying a configuration equivalent to filling the SNO cavity 
with a 1% loaded liquid scintillator or approximately 10 t of isotope after the present 
heavy water experiment is completed. The group is currently researching the optical 
properties of potential nano-crystal solutions. In particular, one must demonstrate that 
sufficient energy resolution is achievable with liquid scintillator. 

6.3.1. Super-NEMO The recent progress of the NEMO-3 program (Sarazin et al 2000) 
has culminated in excellent i3j3(2v) results. In particular, the energy spectra from 100 Mo 
contain nearly 10 5 events and are nearly background free. These data permit, for the 
first time, a precise study of the spectra. In fact, there is hope that the data (Sutton 
2004) will demonstrate whether the I3f3{2v) transition is primarily through a single 
intermediate state or through a number of states (Simkovic et al 2001). The detector 
consists of several thin foils placed between Geiger-drift cells, surrounded by a scintillator 
calorimeter. 

NEMO-3 began operation in February 2003 with several isotopes, 100 Mo being 
the most massive at 7 kg, and plans to operate for 5 years. The collaboration plans 
to increase the mass of 82 Se from 1 kg to 20 kg and begin an additional 5-year run. 
Presently, a Rn-removal trap is being installed to reduce the background, and operation 
should begin again by the summer of 2004. The anticipated sensitivities for TfY 2 are 
5 x 10 24 y and 3 x 10 25 yr for Mo and Se respectively. For the Se data, this corresponds 
to (m pi3 ) below 100-200 meV. 

A much bigger project is currently being planned that would use 100 kg of source. 
The apparatus would have a large footprint however and the Frejus tunnel where NEMO- 
3 is housed would not be large enough to contain it. Currently the collaboration is 
studying the design of such a detector. 
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6.3.2. XMASS The XMASS collaboration (Suzuki 2004; Moriyama 2001) plans to 
build a 10 t natural Xe liquid scintillation detector. They expect an energy resolution 
of 3% at 1 MeV and hope to reach a value for Ty 2 > 3.3 x 10 27 y. This detector would 
also be used for solar-neutrino studies and a search for dark matter. 

6.3.3. Borexino CTF In August of 2002, operations at the Borexino experiment 
resulted in the spill of scintillator. This led to the temporary closure of Hall C in 
the Gran Sasso Laboratory and a significant change in operations at the underground 
laboratory. As a result, efforts to convert the Counting Test Facility (CTF) or Borexino 
itself into a PP(0u) experiment (Bellini 2001; Caccianiga and Giammarchi 2001) have 
been suspended (Giammarchi 2004). 

6.4- The Search for Decays to Excited States 

Searches for /3(3 to excited states in the daughter atom have been performed in a number 
of isotopes but only observed in 100 Mo. (The experimental situation is reviewed by 
Barabash 2000). These experiments typically search for the 7 rays that characterize the 
excited states and therefore are not mode-specific searches. The interpretation therefore 
is that the measured rate (or limit) is for the f3(3(2v) mode. These data may be very 
useful to QRPA nuclear theory because the behavior of the nuclear matrix elements with 
respect to g pp for the excited state decays is different than for transitions to the ground 
state (Griffiths and Vogel 1992, Aunola and Suhonen 1996, Suhonen 1998). Thus, the 
excited state transitions probe different aspects of the theory and may provide insight 
into the physics of the matrix elements. 

A further reason for interest in decays to the excited state, as mentioned earlier, is 
the potential ability to discover the process mediating the decay (Simkovic and Faessler 
2002, Tomoda 2000). However, the decay rate to an excited state is 10-100 times smaller 
than rate to the ground state (Suhonen 2000a, 2000b). Furthermore the structure of 
the excited state in the daughter nucleus is not as well understood as the ground state, 
and this increases the relative uncertainty in the nuclear matrix element. 

6.5. The Search for j3 + j3 + Modes of Decay 

The f3 + modes of decay have not received the attention of the /3~ modes because of 
the greatly reduced phase space and corresponding long half-lives. However, their 
detection would provide additional matrix-element data. Furthermore, if the zero- 
neutrino mode were detected, it might provide a handle on whether the decay is 
predominantly mediated by a light neutrino or by right-handed currents (Hirsch et 
al 1994). 

Radiative neutrinoless double electron capture is a possible alternative to traditional 
neutrinoless double beta decay (Sujkowski and Wycech 2003). In this process, two 
electrons are captured from the atomic electron cloud and a radiated photon carries the 
full Q-value for the decay. A resonance condition can enhance the rate when the energy 
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release is close to the 2P-1S energy difference. In this case, high-Z, low-energy-release 
isotopes are favored (e.g. 112 Sn). Unfortunately the mass differences for the candidate 
isotopes are not known precisely enough to accurately predict the overlap between the 
two energies. If a favorable overlap does exist, however, the sensitivity to {mpp) might 
rival that of /5/9(0z/)decay. 

6.6. Towards a 100-kg experiment 

The KDHK spectrum shows a feature very close to the /3/3(0v) endpoint. This intriguing 
result will need to be confirmed or refuted experimentally. One can see the required 
operation parameters for a confirmation experiment from Klapdor-Kleingrothaus et al 
(2004). One needs about 75 kg-y of exposure, and a background lower than about 0.5 
counts/(kg y). Note that most of the proposals described above will all accomplish this 
very early on in their program if they meet their design goals. If instead one designs 
an experiment only to test the claim (not to provide a precise measurement of the Ty 2 ) 
then a 100-kg experiment could provide the answer after a modest run time. 

If the KDHK result holds up, it will be a very exciting time for neutrino-mass 
research. A {m^p) near 400 meV means that /3-decay experiments and cosmology will 
be sensitive to the mass. As a result, one can certainly imagine a not-too-distant future 
in which we know the neutrino mass and its Majorana-Dirac character. Towards this 
goal, a precision measurement of {mpp) will be required. To accomplish this, we will need 
more than one (3(3 experiment, each with a half-life measurement accurate to 10-20%. 
At this level the uncertainty will be dominated by the matrix element uncertainty even 
if future calculations can be trusted to 50%. With two experiments utilizing different 
isotopes, one might disentangle the uncertainty in |Mo„|. 

6.7. Towards the 100-ton experiment 



The next generation of experiments hopes to be sensitive to JSml tm . If they fail to see 



(3(3(0u) at that level, the target for the succeeding generation of efforts will be JSrn^ 



This scale is an order of magnitude lower and hence will require two orders of magnitude 
more isotopic mass, approximately 100 tons of isotope. 

A 100-ton experiment will have to face the same technical challenges associated 
with radioactive backgrounds and energy resolution as today's proposals. (See Elliott 
and Vogel 2002 for a discussion of these issues.) In addition, a background from solar 
neutrinos will also have to be considered. Solar neutrinos can result in a background 
via elastic scattering or charged current interactions. The rate (Rpp) of (3(3(0u) events 
can be written 



where MW(g) is the molecular weight of the (3(3 isotope and M is the mass of the target 
in tons. For pure 136 Xe, this would result in about 3 /3(3(0v) events per year per ton for 





1 dN _ XN r 420 / 10 2 ' 
M~dt ~W~ MW{g) I ~7j 




(24) 
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T% = 10 27 y. 

Elastic scattering (ES) rates from 8 B solar neutrinos can be comparable to this 
(5(5{Qv) rate if the target material contains the isotope at a low fraction. For example, 
a 2% solution of Xe in liquid scintillator has been discussed as a possible (5(5(0v) 
experiment. Of course the number of ES events within the (5(5(0u) window depends on 
the resolution and therefore we need the cross section per unit energy (Aa/AE)(Bahcall 
1989) 

^«9x 10- 48 cm 2 /keV . (25) 
AE 

The rate of ES events (R 8 ) is then given by 

« (8 x 10-V(keV y t)) (aBjJ^W.) (26) 

where F 8 is the 8 B neutrino flux (5 x 10 6 /(cm 2 s)), Na is Avagadro's number, AE is 
the (5(5(0u) energy window in keV, MW t is the target molecular weight, and N e is the 
number of electrons per molecule of the target. For a pure Xe target with an energy 
window of 50 keV, we find R$ ~ 0.02/(y t). This background is not a problem for any 
pure Xe detector that proposes a half-life sensitivity of 10 28 y. It is significant for a 
detector with only 2% Xe at T^ 2 =10 27 y. 

Charge-current (CC) scattering of solar neutrinos, especially the large flux of solar 
pp neutrinos may also be a background for (5(5{Qv) decay. As pointed out by Raghavan 
(1997), some (5(5 isotopes make interesting targets for pp solar neutrino experiments 
because the reaction produces a radioactive isotope, the intermediate nucleus in the (5(5 
process. The decay of this product nucleus provides a coincidence signature for the CC 
reaction. However for several reasons, this process must be considered as a background 
for (5(5(0u) decay. First, the /3-decay Q-value of the intermediate nucleus is larger than 
the (5(5 Q-value for most of the nuclei. Furthermore, the half-life of the intermediate 
nucleus is often too long for an effective coincidence identification. Finally, the end 
product /3-decay daughter nucleus is the same as the (5(5 daughter. 

The rate of CC events (Rcc) can be written: 

Rcc = Foe Mj^Se = 18 R SNUJ JjL- )( y-' t ->) , (27) 

where Fee is the solar neutrino flux, (ace) is the spectrum-weighted CC cross section, 
f E is the fraction of intermediate nucleus (5 decays that fall within the (5(5(Qv) energy 
window, and Rsnu is the rate of CC interactions in SNU (1CT 36 interactions/ (s target 
atom)). For example, the total R cc ~ 120/(y t) (that is with f E = 1) in the MOON 
proposal is much higher than Rpp m 4/(y t) for a half-life of 10 27 y. In that case however, 
the intermediate nucleus has a convenient lifetime of 16 s, short enough that a delayed 
coincidence can identify the decay and separate it event-by-event from (5(5(0v) candidate 
events. 
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Rcc i n eacn (3(3 isotope must be considered separately. For example, 76 Ge, 116 Cd, 
and 130 Te do not have a low-lying intermediate nucleus level reachable by the high-flux 
pp neutrinos, and in the case of Ge, 7 Be neutrinos. Alternatively, 82 Se, 176 Yb, and 
160 Gd have significant cross sections but the long-lived intermediate nucleus prevents 
easy anti-coincidence identification. In the important case of 136 Xe, little is known 
about the structure of the intermediate-nucleus ( 136 Cs) excited states. The decay of 
136 Cs is to highly excited levels in 136 Ba that decay themselves via 7 emission. This 
results in a relatively large but also a 7-ray cascade that might provide a signature 
to eliminate the potential background. 150 Nd is another isotope in which little is known 
about the intermediate-nucleus excited states. 

7. Conclusion 

If the reader retains anything from this review, it should be that the information 
recently acquired from oscillation experiments makes this an exciting time for (3(3 decay 
because th e next generation of experiments will be sensitive to neutrino masses on the 
order of \J $rn? &tra . If a nonzero rate is seen, we will know that neutrinos are Majorana 
particles. With some progress on the calculation of nuclear matrix elements — and 
we believe progress is possible — a nonzero rate in these experiments should allow the 
determination of the hierarchy realized by nature and the absolute mass scale. If no 
signal is seen, we should be able to say either that the hierarchy is normal or that 
neutrinos are Dirac particles. 

Although none of the next-generation experiments is ready to operate, many of the 
new proposals are promising. We are hopeful that with concentrated effort from theorists 
and experimentalists, (3(3 decay will add to our growing understanding of neutrinos. 
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